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When fast cracks become unstable to microscopic branching (micro-branching), fracture no longer 
occurs in an effective 2D medium. We follow in-plane crack front dynamics via real-time measure¬ 
ments in brittle gels as micro-branching unfolds and progresses. We first show that spatially loeal 
energy balance quantitatively describes crack dynamics, even when translational invariance is badly 
broken. Furthermore, our results explain micro-branch dynamics; why micro-branches form along 
spatially localized chains and how finite-time formation of cusps along the crack front leads to their 
death. 
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Although fracture mechanics is a mature field, many of its simplest questions remain unanswered. Cracks govern 
material stability. The classic 2D descriptions of fracture [1] shows that the existence of a crack focuses elastic energy 
stored in a strained material into a single point - the crack tip. There, stresses diverge as a ^ K|^/r where r is the 
distance from the tip and AT is a coefficient called the stress intensity factor. The motion of the crack is then governed 
by the physical demand of energy balance; that G, the elastic energy flowing to the crack tip, equals the fracture 
energy T (the energy dissipated per unit area of crack extension). For example, a sudden increase in T will cause a 
corresponding decrease in crack velocity, since part of the elastic energy that is used to separate the crack faces is 
now used to compensate for the increase in fracture energy. Energy balance is the basis for the continuum description 
of the dynamics of simple, straight cracks in brittle materials. 

Cracks, however, fracture 3D materials. Whereas a crack’s tip in 2D materials is a singular pointy in 3D the leading 
edge of a crack forms a singular /me, the crack front In “simple” cracks, the crack front is a straight line with no overt 
dynamics. For this reason, 2D descriptions of fracture dynamics for simple cracks are sufficient and the classical 2D 
theory delineated above [1] describes crack dynamics perfectly [2-4], as long as they propagate along a straight-line 
trajectory. 

Simple cracks are generically unstable entities [5] that leave complex 3D surfaces in their wake [6-9]. Their intrinsic 
3D character calls for a study of the dynamics of crack fronts [10-14]. Even when a crack is constrained to a plane, 
the geometry of the crack front is important. For example, the advance of quasi-static fronts in patterned materials 
is influenced by the long-range elastic forces that act along crack fronts [15-17]. A recent study [18] of quasi-static 
peeling of adhesive tape demonstrates the striking consequences of manipulating front geometry to enhance material 
toughness. The elasticity of crack fronts is also important to explain the fluctuating dynamics of quasi-static planar 
cracks propagating through random heterogeneities [19-21]. 

When cracks are rapid, experiments suggest [22] that distortions of a crack front produce inertia in the “massless” 
cracks described by 2D theory; an intrinsic 3D effect. The micro-branching instability of rapid simple cracks is 
an example where a crack front loses its symmetry dynamically. Simple cracks in brittle materials experience this 
instability above a critical crack velocity Vc ^ 0.3ci^, where cr is the material’s Rayleigh wave speed [5, 22, 23]. 
For velocities v > Uc, directed chains branch-lines^^) of micron-scale branched cracks are spontaneously generated 
(see Fig. lb). Each “micro-branch” branches away from the main crack to propagate beneath its fracture surface. 
Every micro-branch creates additional fracture surface and hence, increases the value of T that is felt by the main 
crack [24]. As, for short times, the total energy flux to the crack front is constant, micro-branching also reduces 
the energy flux to the rest of the crack. A single micro-branch may effectively increase T by up to 100% (if it 
creates an additional crack propagating parallel to the main crack). One may therefore consider micro-branches as 
energy sinks that are dynamically “toggled” on and off. As micro-branches are localized along the crack front, these 
perturbations are spatially localized in both parallel and transverse to the propagation direction. Therefore, when 
excited, micro-branches locally perturb initially straight crack fronts to produce large fluctuations in v. 

In this Letter, we study the propagation of fast crack fronts during the micro-branching instability. We will provide 
a detailed account of in-plane front dynamics as micro-branches nucleate, grow and eventually die. This will provide 
us with new insights as to how front geometry and motion continuously feed each other; eventually conspiring to 
cause micro-branch “death”. 

We study crack front dynamics by performing real-time visualization of the fronts as they traverse the fracture 
surface. We do this by using brittle polyacrylamide gels composed of a 13.8% (w/v) of acrylamide and a 1:37.5 
(w/w) bisacrylamide-to-acrylamide ratio where cr = 5.2 m/s. When scaled by cr, these gels exhibit both the same 
single crack dynamics [4] and micro-branching phenomenology [23] as more conventional brittle materials such as 
glass and PMMA. The advantage of using gels is in reducing the wave speeds, hence crack velocities, by three orders 
of magnitude. 

Our experimental system is schematically described in Fig. 1(a). Our gel samples are cast to be rectangular blocks 
of dimensions 50 x 60 x 6mm^ {X x Y x Z), where X, T, and Z are, respectively, the propagation, loading and 
thickness directions. We image the crack front by shining collimated LED light through the samples in a direction 
normal to the XZ (fracture) plane. The strong curvature at the crack tip deflects the light at the crack front, creating 
a sharp shadow boundary in the image plane and producing a 2D projection of the crack front. We capture front 
dynamics by imaging the front using high-speed (IDT-Y4) camera at ^48000 frames/sec in a “X x Z” window of 
10 X 1 mm^ located at the center section of the gel sample. Our spatial resolution was ^ 10 micron per pixel. 

We generated crack front dynamics in the micro-branching regime by displacing the system’s boundaries by a 10- 
12% strain, prior to initiating fracture. Fracture was initiated by inserting a small “seed” crack at the sample’s edge, 
midway between its vertical boundaries. This procedure produced crack velocities of 0.1-0.5ci^ along the mid-plane 
of the sample. Ah velocities, v(z), referred to in this paper are the normal velocities to the front at each point z. 

In Fig. 1(b) we present a typical picture of the fracture surface created. Crack propagation is complex, forming 
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FIG. 1. (a) The experimental setup. A rectangular block of polyacrylamide gel is strained by displacement of its Y boundaries. 
Crack fronts, propagating within the sample’s mid-plane are visualized via shadowgraph with a high-speed camera; illumination 
is via collimated LED light normal to the fracture (XZ) plane, (b) A post-mortem fracture surface section containing four 
nearly parallel directed chains of micro-branches (“branch-lines”), (c) Top. The imaged fronts that formed the branch-line 
denoted by the box in (b). Bottom. A close up photograph of the corresponding fracture surface with overlaid edge-detected 
fronts. 


micro-branches at several locations along the front. The four branch-lines appearing in the figure are oriented parallel 
to the mean crack velocity. To characterize the corresponding front dynamics, we focus on the formation of a single 
branch-line, denoted by the dashed box in Fig. 1(b). In Fig. 1(c) we present a series of snapshots of successive 
crack fronts within this box, together with a close up of the resulting post-mortem branch-line formed by these fronts. 
A direct comparison between the fronts and resulting fracture surface (Fig. 1(c), bottom) shows that a one-to-one 
correspondence exists between micro-branches on the fracture surface and loci of high curvature along the front. 

Fig. 2 presents a detailed picture of the motion within a typical section of a front that includes three well-defined 
micro-branching events. Along the crest of a single branch-line marked A we can see that, prior to the nucleation of 
the first micro-branch, ( 1 ) the front has nearly zero curvature. Then, with the nucleation of the first micro-branch 
(2) the local front velocity along A gradually decreases. The apparent cause for the slowing down is the effective 
increase in local fracture energy due to the incipient micro-branching. The resulting gradient in velocity along the 
front produces an increase of front curvature. The front acquires a locally concave shape which becomes more and 
more pronounced until (3) the curvature suddenly drops as the local velocity peaks to ^ 100% above the mean 
front velocity. Examination of the fracture surface suggests that the release of curvature and the velocity jump are 
coincident with the “death” of the micro-branch. Immediately afterwards (4) the velocity drops again as a new 
micro-branch nucleates. In contrast to the strong fluctuations along A, the “quiet” adjacent mirror-like region [25] 
marked B shows practically no change in curvature. We do, however, see significant variations in velocity which, after 
a delay, follow the changes in velocity along the branch-line A . 

Let us consider these results. First, the association of micro-branch nucleation and death with local velocity decrease 
and increase along the branch-line is qualitatively consistent with local energy balance; G = F. Our second observation 
is that after the front acquires curvature and the micro-branch dies, the local velocity does not only return to the 
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FIG. 2. Three consecutive branching events along a branch-line, (a) The series of fronts captured during branch-line formation. 
Colors indicate local velocity levels. Marked points show ( 1 ) an unperturbed front ( 2 ) local velocity decrease associated with 
micro-branch nucleation (3) maximum local curvature (4) formation of a new micro-branch, (b) Plots of local curvature ( Top) 
and local velocity (Bottom) along the cuts A (dots) and B (open symbols) in (a), (c) The corresponding post-mortem showing 
the branch-line associated with the fronts in (a). 


mean velocity level, but doubles it. These observations suggest that front curvature must be increasing G locally to 
produce these high velocity peaks. This idea is supported by a well-known result derived by Rice [26] for static fronts. 
Rice computed the static stress contribution made by a small in-plane perturbation to a straight front. When a front 
X = xq -\- Sx(z) is perturbed in the sense that \dx/dz\ <C 1, the stress maintains the inverse-square-root singularity at 
each point along the front, but K = Kq T- 6K(z) varies along the front as 

Ko 2 tt ^ J (z' - ^)2 ’ ^ ^ 

where xz is the fracture plane, x is the crack propagation direction and the unperturbed front is x(z) = xq with 
a stress intensity factor Kq. A salient feature of this integral is that it is solely determined by the front geometry. 
It generates a stabilizing force (“line tension”) that tends to restore the front to a flat configuration, since G oc 
or 6G/Go = 26 K/Kq. This force is countered by the local increase of the fracture energy, (^r( 2 ;), generated by the 
micro-branch. In general, we cannot measure ^F directly. When, however, a micro-branch “dies”, (5F ^ 0 and the 
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elastic energy stored in the front curvature locally accelerates the crack. Evidence for this is seen Fig. 2b along the 
A line, where we see that each peak in curvature is followed by a peak in local front velocity as the line tension is 
released. 


a 


2bK(z)IK. 6v(z) [m/sl ^ 1 .2 



FIG. 3. Comparison of the instantaneous velocity and local stress intensity factor, computed by Eq. (1), along a front at the 
moment of micro-branch “death”. (a) Spatial stress (dots) and velocity variations (open circles) for the front in the inset are 
well-correlated, {inset) Overlay of the front at the moment of maximal velocity on the resulting post-mortem fracture surface. 
Front positioning is approximate, (b) Sv{z) vs SGjGo = 26K/Kq for both the data in (a) and twelve such releases from other 
events. Different colors indicate different events. The dashed “y=x” line is a guide to the eye. 


While these statements appear qualitatively accurate, are they quantitatively correct? Let us now focus on the 
moment of micro-branch death; the moments corresponding to the peak in local front velocity at the instant where 
(5F ^ 0. We consider thirteen cases, with mean crack velocities ranging between 0.1 — 0.3ci^, where we identify a clear 
transition from a branched to smooth surface in the post-mortem image. 

At these instances we can directly compare the relative velocity deviation from the mean 6v{z) = v(z) — {v{z)) and 
the built-up energy release rate 6G/Go = 26K{z)lKo computed by Eq. (1). As Eq. (1) was derived for fronts parallel 
to the 2 : axis, we computed 6K/Kq over symmetric regions of the front around the point of maximum curvature by 
first rotating the front to make edges of this region sit on an x = const line. The average {v{z)) was computed for 
velocity values in the same region. 

The result of the comparison is presented in Fig. 3. We see that the correction for the stress intensity factor given by 
Eq. (1) is clearly correlated with the instantaneous velocity profile at the moment of micro-branch “death”. Despite 
being only a first-order correction, Eq. (1) is shown here to work quite well. This is demonstrated locally along a 
typical front in Fig. 3a and in another twelve similar instances in Fig. 3b. It is a manifest demonstration of how front 
curvature is translated into motion. We note that we see no observable effect of the proximity of neighboring branch 
lines in this collapse. This is consistent with the observations of [27] in quasi-static propagation. Since we compare 
here a dimensionless quantity with velocity Sv = 2vo6KIKq, the coefficient of proportionality vq must also have the 
dimensions of velocity. Its values vq = 1 ±0.3m/5 are far from the characteristic wave velocity {cr ^ bm/s), but close 
to the measured mean crack velocities {v) = 0.7 — l.bm/s. This might be related to the strong velocity dependence 
of the fracture energy in polyacrylamide gels [4]. We give a possible origin of this coefficient in the supplementary 
material [28, 29]. 

Let us now consider the front evolution depicted in Fig. 4 where the process leading to the build-up of the “line 
tension” prior to the death of a micro-branch is highlighted. During the buildup of the line tension, v along the crest of 
the branch-line is approximately constant during the rapid increase in curvature that leads to cusp formation. In Fig. 
4(a) we present a series of fronts leading to a release event. The release is preceded by a rapid increase in curvature 
which culminates when the front attains a cusp-like form. As in the events depicted in Fig. 2, local increases in v 
immediately follow cusp formation. 

It is well-known that propagating fronts may develop shocks in curvature [30], or cusps. If a cusp does form, we 
would expect the stress intensity factor to locally diverge according to Eq. (1). It is, therefore, likely that this high 
stress could promote the “death” of the micro-branch and hence trigger the local release of a front. 
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FIG. 4. (a) Cusp formation during a typical micro-branching event. We present every second front. Red arrows denote normal 
displacements, (b) Time-series of the peak value bKjK^ (fop) and front velocity (bottom) along the crest of a branch-line show 
a few cycles of build-up and release. We dehne the micro-branch lifetime, r, by the interval over which SKjKo hrst surpasses 
(micro-branch birth) and then drops below (micro-branch death) a threshold value 0.3 (red line). Note that during stress 
build-up, the local velocity stays approximately constant. Shaded areas: 3 examples of r for different branching events, (c) r 
compared to predicted cusp formation times = l/{KoVn)- ko is the initial maximal curvature and Vn is the initial normal 
velocity at the location of maximum curvature. Colored points denote events in (b) and those shown in the inset, {inset) 
Curvature blow-up dynamics. Time is normalized by t*. t = 0 corresponds to the initial front, n is the maximal curvature of 
each successive front normalized by kq. Black dashed line: predicted finite-time blow-up. 


How can this picture be tested? A nearly constant normal velocity, v{z^t) = where Vn is taken to be the initial 
velocity at the beginning of the event, coupled with an initially curved section of the front will spontaneously produce a 
cusp with no further assumptions. To see this, we rewrite v = Vn terms of front coordinates; dtx = {dzx)‘^. 

Differentiating by we obtain an equation of motion for the local slope u = —dzX: 


du 


du 




= 0 . 


( 2 ) 


It is known that this equation develops shocks (i.e. cusps in u) in finite-time where d^x = —dzU ^ (C — and 
= l/{f<iQVn) when /^o is the maximal front curvature at t = 0. 

The constant velocity model provides us with a prediction for the time required for the front to form a cusp C. 
Is this timescale related to the lifetime, r, of a micro-branching event? To test this, we consider an ensemble of 
micro-branching events. As depicted in Fig. 4(b), we define r as the time between the initial stress build-up and 
stress release points for each event, using a threshold of dKjKQ = 0.3 to define these points. We use the initial 
maximal curvature, /^o, and the normal velocity, Vn^ to evaluate the predicted cusp formation time, t*. (We measure 
hq by rotating the front to eliminate any mean slope and fitting the region of highest curvature with a parabola, 
is the average normal velocity over the same range in 2 :.) We consider only events where the estimated error in C was 
no larger than 50%. 

As shown in Fig. 4(c) this simplistic model works rather well. Within the limits of our experimental resolution 
T over the stress build-up we observe rapid growth of curvature that culminates near the predicted Moreover, 
our predicted C indeed correspond to the death of each micro-branch. 
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This being said, micro-branch dynamics do not always proceed directly to curvature blow-up. In some cases the 
front may “slip”; undergoing partial release during an event before final cusp formation. An example of such slip is 
the curve denoted by open symbols in the inset of Fig. 4 (c). 

In conclusion, we have presented the first experimental exploration of rapid crack front in-plane dynamics. Our 
observation that micro-branches act as sinks that effectively increase the fracture energy agrees qualitatively with 
a local energy balance. The induced front curvature may explain the micro-branch localization in z as the velocity 
overshoots upon micro-branch release appear to generate new micro-branches, thereby forming branch-lines. Moreover, 
we have shown that the pronounced velocity overshoot at the moment of micro-branch “death” is well correlated with 
the static contribution for the stress intensity factor resulting from front curvature. This is a clear demonstration of 
how stresses distributed along a crack front translate locally into dynamics. 

Nevertheless, the demonstrated validity of the expression (Eq. (1)) is intriguing for several reasons. First, Eq. (1) 
is a first-order correction in the local slope of the front, while the fronts we considered contained mild slopes. Our 
current data do not have the precision to investigate the effects of higher order corrections [31-33] and it would be 
interesting to test their implications in the future. Secondly, Eq. (1) was derived for static fronts while we consider 
moving fronts. In the work of Ramanathan & Fisher [10] the static theory was expanded to include dynamic effects. 
These include waves that propagate on the front at velocity ^ — v‘^. Since in our experiments v‘^ ^ 0.1 ^ 1 

and because of the small size of micro-branches ( ^ 100//) and relatively long life-time IOO// 5 ) we expect any 
signature of the waves to average during micro-branch growth. We might expect the agreement with Eq. (1) to break 
down at larger vjcR, when such inertial effects become important. Lastly, apart from the increase in fracture area 
caused by the out-of-plane motions of the front, micro-branches should also influence the local energy release rate by 
introducing shearing stresses. The combined effects of these contributions on dynamics demand further study. 

In addition, while recent experiments [34] provide a mechanism for micro-branch birth, we have seen that the front 
dynamics induced by this “birth” could lead to eventual micro-branch “death”; front curvature caused by micro¬ 
branch initiation may lead to micro-branch death via the large line tension generated by dynamic cusp formation. 
This scenario could provide a key to understanding why micro-branches remain small and do not develop. The above 
results are all purely 3D effects that underline the necessity of extending fracture mechanics to the third dimension. 
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fronts. 
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Supplemental Material for Crack Front Dynamics 


A HEURISTIC CALCULATION OF vq 


From Fig. 3 of the main text we concluded that at the moment of micro-branch “death”, the local velocity is 
proportional to the line tension given by Eq. (1). The origin of the proportionality constant, needs still to be 
accounted for. Here we suggest a possible heuristic explanation. 

Our statring point is the energy balance between the energy release rate and the fracture energy 

G = V. (SI) 


We proceed by locally perturbing both sides of the balance to take into account the variation of F and G with local 
front geometry and velocity. First, we note that G oc g{v)K‘^ where g{v) l—vjcR is the dynamical correction factor 
of G [SI] and K here is the non-dynamic part of the stress intensity factor. Then perturbing around a straight front 
configuration with a constant energy release rate G = Go we have 


^G 


d\ogg 

dv 


6v{z) + 2 


6K{z) 

Ko 


(S2) 


where 6K{z)/Ko is given by Eq. (1). Next, let’s discuss the fracture energy. In polyacrylamide gels the fracture 
energy for a simple crack is a strong function of velocity [S2]. We will term it here To{v). Since in the instances 
we consider in Eig.3 there is no contribution from the micro-branch to the fracture energy, the local perturbation is 
simply 


^F dlogFo 
i 0 dv 

Equating SG/Gq = ^F/Fq we find that Sv = 2vo6KlKo with 

'dlogFo dlog^ 

_ dv dv 

To evaluate the expression for vq we use the results of [S2, S3]. We approximate the fracture energy by a linear function 
ro('i;) ^ T{vm)-^^o{vm){v-Vm) arouud Vm = O.Scr. Then d\ogTo/dv = l/{v^u) where u = V ^{vm)- Vm = 
(0.25 =b 0 .1)ci^. Noting that also d\ogg/dv zz —1 /[cr — v) we see that vq = {v u){l — v/cr)/{1 -h u/cr). Since the 
mean front velocities for the data presented in Fig. 3 vary in the range v = 0.8 — 1.6m/5 the calculated proportionality 
coefficient is then vq = 1.4 ± 0.3m/5, which is within the margin of error of the value of vq obtained from Fig. 3. 


(S4) 


VALIDITY OF THE CONSTANT VELOCITY ASSUMPTION DURING MICRO-BRANCH GROWTH 

In the main text we argued that during the stress build-up phase of a micro-branching event the velocity stays 
virtually constant in the region of the micro-branch. Here in Fig. SI we present statistics of 37 events where 
we compute the normalized standard deviation of local velocities for each events. Events where velocity markedly 
deviates from the mean are rare. 
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FIG. SI. The constant velocity assumption for the cusp formation model, (a) An example of a series of consecutive crack 
fronts. In red are the ’’cusp” regions where the velocity statistics were taken from, (b) A histogram of normalized standard 
deviation values for the 37 micro-branching events. For each event considered in Fig.4c of the main text, we gathered velocity 
values from all the ’’cusp” regions of all the fronts, taking the fronts from the first passing the threshold to the one reaching 
maximum stress. Then for each event we calculated the mean value of all velocity values {v)event and used it to normalize the 
velocity values. 




































































